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Abstract 

We measured the current induced loss of vortex correlation in YBCO crystals 
using the pseudo DC flux transformer geometry. We find that the current Icut 
at which the top and bottom voltage drops differ is a linear function of tem- 
perature independent of pinning. Although the vortex correlation length at 
Icut coincides with the sample thickness, the experimental data show that I^ut 
is sample thickness independent. The observed behavior with temperature, 
magnetic field and thickness is described by a phenomenological model. 
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The laminar structure of the high temperature superconductors has important influence 
on the vortex behavior in the mixed state. This laminar structure and the presence of ther- 
mal fluctuations yields an extremely rich H — T phase diagram.lll Recent experiments!, using 
a DC flux transformer contact configuration, have shown that the correlated vortex motion 
can be destroyed by either thermal fluctuations or the application of large driving forces. 
The non-homogeneous current distribution induced in electrical transport measurements 
in YBCO(123) single crystals, in the DC-transformeiS configuration, has been usedHi'i to 
determine the correlation length of the vortices in the direction of the applied field (c crys- 
tallographic direction). It was shownii that there is a well defined temperature, Tth{H), 
above the irreversibility line Ti{H), where the correlation of the vortex velocity across the 
sample is lost. This has been interprete dii decoupling of the vortices nucleated in 
the Cu-0 planes. Since this decoupling takes place in the linear current-voltage response 
regime it is concludedi'i that at that temperature the vortex correlation length in the di- 
rection of the applied field coincides with the thickness, d, of the sample. As a consequence 
the transport measurements in samples of different thicknesses have been usecfl to deter- 
mine the temperature dependence of the correlation length in the c direction at different 
fields. On the other hand transport measurements in the linear regime using a Montgomery 
type analysis have led to the conclusion! that the transport properties of the vortex liquid 
phase has a non-local character.l The non-local character of the electrodynamics induces 
an enhancement of the apparent anisotropy of the resistance of the vortex liquid in the DC 
transformer configuration. ^ 

In this paper we concentrate our interest in the study of the I-V characteristics for T < Tth 
and temperatures above and below the reversibility line, defined as the temperature Ti{H) 
where the resistance of the ah planes tends to zero. For T > Ti{H) the I-V characteristics 
showi a linear regime at low currents and become non linear at higher currents. For T < 
Ti[H) the response is non lineaii for all currents. It is shown that for T < Tth flux cutting 
is induced by a current Jcut{H,T,d), in the non-linear response regime. 

The experimental results are consistent with a current distribution flowing; near the sam- 
ple surface. This allowed us to extend the model developed by J. Clem et alB to explain our 
results obtained using a non-homogeneous current distribution. The results of the model 
provide a qualitative comprehension of the field, temperature and sample thickness depen- 
dence of JcutiH,T, d). It also shows that the understanding of the thickness dependenceSi 
of Tth goes beyond the assumptions made by the model and requires the analysis of the 
contribution of thermal fluctuations in the vortex transport properties. 

The YBCO single crystal was prepared as indicated in Ref. 5. The results presented 
in this paper were obtained using the same samples and contact configuration described in 
Refs. 5 and 6. 

In Fig. 1 we plot typical I-V characteristics for temperatures above (la) and below (lb) 
Tj = 89. 77^ at a field of lOkOe. In this case the sample thickness is d = 20/xm. The 
experiment was made injecting current through contacts 1 and 4, /14, and measuring the 
voltages V23 and Vgy, see insert in Fig. 1. At low currents V23 = Vqi indicating that the 
vortices move correlated across the sample. Increasing the driving current, V23 becomes 
different from Vq7 at a well defined current Icut (arrows in Fig. 1) where flux cutting is 
induced. 

Figure 2a shows the cutting current density defined as Jcut = Icut/{d x w) (where w is 
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the width of the sample) as a function of temperature for samples of different thickness. 
The data show that J cut decrease linearly with temperature vanishing at T^/i.iil For a given 
temperature, Jcut decreases when increasing the thickness. In Fig. 2b we show the same 
data but scaled by the sample thickness. It is clearly seen that Jcut x d becomes thickness 
independent. This result indicates that flux cutting is induced by a force that can be 
associated with a surface current density given by Icut/w. That is, the transport current 
should be considered as flowing in a region close to the surface in a thickness much smaller 
than d, in agreement with the suggestion made in ref. 7 and with the experimental results 
obtained in the linear response regime.! 

It is interesting to point out that the data above and below Ti{H) fall on the same 
straight line. This is demonstrated in Fig. 2a where Tj {Ti/Tth = 0.985) is indicated by the 
arrow for the sample with d = 48/im. Thus, flux cutting should be related to the properties 
of the flux structure, independently of the presence of effective pinning centers. That is, the 
dynamic vortex correlation length in the field direction is not modified by the introduction 
of pinning. 

With the results and considerations made before we present now a simple model which 
enables to rationalize the observed field and temperature dependence of Icut- The model 
also explains why Icut is independent of pinning whereas the voltage Vcut at the cutting does 
depend on pinning. 

Consider a stack of N superconducting layers with vortex pancakes moving in these 
layers. The equation of motion for vortex number i in layer number u (see Fig. 3a) is given 



assume that the effect of the external current is to give rise to a driving force /, which is 
confined to the top layer. This assumption is consistent with the experimental observation as 
discussed above. It is also in accordance with the phenomenological viscous hydrodynamics 
of Huse and Majumdar.l However, it should be mentioned that Huse and Majumdar work in 
the limit of linear resistance, whereas the cutting by current always occur in the non-linear 
region. 

We want to determine the average velocity in layer u: Vy = (fi,^), since the voltage drop 
over that layer is proportional to v,^. We consider the situation where the average velocity is 

the same in each layer, i.e., v = vi,V2, vn, thus Vtop = Vbottom(y23 = Vqj). Assuming that 

the average of the pinning force over vortices and time is independent of the layer number, 
we obtain the following equation: 



by: 




(1) 





(2) 



This set of equations is readily solved to give: 
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The largest stress is between the two top layers. The driving force at which the velocities 
in layer 1 and 2 starts to differ is given by: 

fcut = ^ (4) 

where f^"-^ is the maximum force which the bonds between pancakes in adjacent layers can 
sustain. Notice that fcut is independent of the pinning forces. In contrast, the velocity (or 
voltage) at cutting does depend on the pinning. By substituting Eq. 4 into the equation for 
the average velocity one obtain Vcut = f"^°-^/{N — 1) — Thus if one knows the friction 
coefficient rj one can determine the intrinsic property f"^"-^ together with the extrinsic pinning 
force by a combined measurement of the voltage and the current at cutting. 

To estimate f"^"-^ is like to estimate the yield strength of a material. It is not easy. Let 
^Xmax denote the displacement of vortices in layer 2 at which plastic slip sets in (see Fig. 
3b and 3c). Let AE denote the associated increase in the energy due to the deformation. 
We estimate 

AE 

rmax ^j-^ , » 



^max 



Slip will occur when Axmax reaches some fraction of Oo, the distance between the vortices. 
For small (Fig. 3b), i. e. large magnetic field, we expect AE = kiAx"^ where ki will 
be a short wave length tilt modulus. For lower fields (Fig. 3c), increases and so do 
Axmax- In this case we should rather connect AE with the energy of the Josephson vortex 
running in between the planes connecting the vortices in the two layers. This suggest that 
AE becomes linear in Ax for low fields: AE = K2AX. Both constants ki and k,2 will scale 
with the superfiuid density The maximum force will behave like 

fmax ^ m'G{B), (6) 

where G{B) = 1 for low fields and G{B) ~ B^^^"^ for high fields. 

We have \'ip\'^ ~ (1 — ;^)(1 — -g^^).0 Furthermore, fmax = (poJcut^y where / denotes 
the thickness of the layer {I <^ d) that carries the transport current which gives rise to the 
Lorentz force. In the model Jcut = Icut/{^ x w), therefore 

hut^{l-^){l-j^)G{B\ (7) 

J-c -Dc2U ) 

independent of /. 

Using Bc2{T) = (1 - J)Sc2(0) we finally obtain 

According to this expression, I^ut depends linearly on the temperature. As B is increased 
the intersection of Icut{T) with the x-axis decreases linearly in B. The slope of Icut{T) is 
independent of B for small B. For large B the slope decreases as B^^^'^ . This is in qualitative 
agreement with the experimental results shown in the inset of Fig. 2b. For a quantitative 
analysis of the results, more data at higher fields are necessary. 
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As was discussed above, the experimental results shown in Fig. 2a prove that Icut is 
independent of pinning. This agrees with our model, as shown by expression (1), where fcut 
is found to be pinning independent despite pinning is included in the equations of motion 
(see equation |ID. Thus, expression (j^) describes the experimental temperature and field 
dependence of Icut and shows that Icut is an intrinsic property of the flux structure, in 
agreement with the data of Fig. 2a. 

We have restricted our discussion to purely mean field considerations neglecting thermally 
induced vortex fluctuations. We now want to turn our attention to the field dependence 
of Tfh. As shown in Fig. 4 one finds experimentally that T^h depends linearly on the 
magnetic field. In accordance with the definition of Tth, it appears natural to associate this 
temperature with the temperature where vortex-antivortex pairs thermally activated in the 
superconducting planes unbind. The unbinding in zero field happens when the ratio T/lipl"^ 
achieves a certain invariant value0: 

Tub/\i'{Tub)\'^ = const. (9) 

It has previously been sugg ested0 that the effect of the external magnetic field on the 
unbinding transition in a first approximation can be obtained simply by including the mean 
field dependence on 5 in l-^/^p in Eq. (y). 
This leads to 

= 1 E . (10) 

The linear field dependence of Tub and the decoupling of the layerslli at T^b, strongly 
suggest that T^t is identical to Tth- The dependence of T^h on thickness arises then as an 
effect of the correlation length (in the field direction) connected with the unbinding transition 
becoming equal to the sample thickness.il 

We have been able to show that flux cutting in YBCO(123), using the DC transformer 
configuration, is induced by currents flowing at the sample surface, in strong agreement 
with the suggestion made in ref. 7 and the previous results of ref. 2, obtained in the linear 
response regime. This has allowed us to formulate a simple model which disregarding the 
effect of thermal fluctuations, describes qualitatively the results in the non-linear response 
regime. The temperature dependence is correctly described by the model, assuming that 
the coupling between planes has its origin in a Josephson energy. 

We thank E. Jagla, C. Balseiro and E. Osquiguil for discussions in the course of preparing 
this manuscript. This work was supported by The British Council - Fundacion Antorchas 
and by the British EPSRC. 
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FIGURES 



FIG. 1. V23 and Vgy as a function of driving current /14 for two different temperatures: a) 
89. 9K and b)89.4K. Also shown is AV = V23 — Vqj. The sample geometry is shown in the inset. 
The arrows show the current, Icut, where V23 and Vqj start to differs. 

FIG. 2. a) Temperature dependence of Jcut = Icut/{d x w) for different sample thickness, for 
an applied field of lOkOe. The arrow shows the corresponding Tj for the sample with d = 48/um. 
The solid lines are linear fit to the data, b) Same data of Fig. 2a scaled by the sample thickness d: 
{Jcut X d) as a function of temperature. In the inset is shown the temperature and field dependence 
of Icut for the sample with d = 32/im. The data were fitted using expression (^) (solid lines). 

FIG. 3. Schematic diagram showing the velocity of the pancakes vortices in different supercon- 
ducting layers (see text). Sketch of the vortex configuration in adjacent layers for large (b) and 
small (c) fields. Ax is the relative displacement between vortices in adjacent layers due to the 
external current. 

FIG. 4. Phase diagram showing the temperature Tj and Tf^ for a sample with d = 20/im. The 
solid line is a linear fit and the dotted one is a guide to the eye. 
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